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. Abstract. The amount and nature of the evolution of the X-ray prop- 

erties of clusters of galaxies provides information on the formation of 
structure in the universe and on the properties of the universe itself. The 
cluster luminosity - temperature relation does not evolve strongly, sug- 
gesting that the hot X-ray gas had a more complicated thermodynamic 
\^ \ history than simply collapsing into the cluster potential well. Cluster X- 

ray luminosities do evolve. The dependence of this evolution on redshift 
and luminosity is characterized using two large high redshift samples. 
QQ ' Cluster X-ray temperatures also evolve. This evolution constrains the 

i dark matter and dark energy content of the universe as well as other 

parameters of cosmological interest. 

o 

O 

1. Introduction 

Q , The evolution of cosmic structure is strongly dependent on the cosmology of the 

\-{ ' universe. Jenkins et al. (1998) is one of the many papers describing this well 

, known result. Clusters of galaxies, as the most massive bound objects known, 

^ ' are the ultimate manifestations of cosmic structure building. The evolution of 

clusters is simple, being driven by the gravity of the underlying mass field of 
the universe and of a collisionless collapse of cluster dark matter. It should be 
possible to calculate this evolution reliably compared to that of other objects 
I visible at cosmological distances such as galaxies or AGN. Clusters are luminous 

X-ray sources. The X-ray emission mechanism is optically thin thermal radiation 
from a medium nearly in collisional equilibrium, about the simpliest situation 
imaginable. Thus observations of the X-ray evolution of clusters provide a robust 
measure of the evolution of cosmic structure and thereby constrain the cosmology 
of the universe. 

The cosmological model is described by a set of cosmological parameters. 
The present value of the Hubble parameter is Hq = 100 h km Mpc~^. When 
needed h = 0.5 will be used, but almost nothing about evolution depends on the 
precise value of h since data at two epochs are always compared. The present 
matter and dark energy densities in terms of the critical density are and 
JIao respectively. The amount of structure in the universe is described by fig, 
the present rms matter fluctuations in spheres of 8 h~^ Mpc. This parameter 
is a complicated way to describe the present normalization of the spatial power 
spectrum of matter density fluctuations, P(k), on a scale of k ~ 0.2hMpc~^: 
crs ~ [P(0.172hMpc-i)/3879h~3]Yipc3]i/2 (Peacock, 1999 equations 16.13 and 
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Figure 1. The cluster luminosity - temperature relation. Triangles 
are the low redshift sample and squares arc the high rcdshift sample. 
The line is the best fit to the entire sample of L{hol) = C[kT]" {1 + z)^ . 
a. (left) For a cosmology with = 1-0, ^ko = 0.0. b. (right) For a 
cosmology with = 0-3, f^Ao = 0.7. 



16.132). The dark energy equation of state is P = w p c^. If w = -1, then the dark 
energy is the cosmological constant, if — 1 < w < it is termed Quintessence. 
Recall that w = is cold dark matter and w = 1/3 is radiation. Cluster temper- 
ature evolution provides constraints on all of the above cosmological parameters 
except h. 

When needed, two specific cases will be considered. The X-ray astronomer's 
universe where fimO = 1-0, J^ao = 0.0. This combination was known to be correct 
ten years ago. The other case may be called the bandwagon universe where 
^mO = 0.3, r^Ao = 0.7. This combination is known to be correct today. Both 
universes are spatially fiat, + ^Ao = 1- 



2. Cluster Luminosity - Temperature Relation Evolution 

The relation between cluster X-ray luminosities and temperatures (L - T rela- 
tion) is the oldest and best studied among cluster X-ray observables. Early stud- 
ies include Mitchell, Ives k, Culhane (1977), Mushotzky et al. (1978) and Henry 
& Tucker (1979). Recent comprehensive work includes David et al. (1993), 
White, Jones & Forman (1997), Markevitch (1998) and Arnaud & Evrard (1999). 
Figure 1 shows an example of this relation using a sample of 53 z < 0.3 and 32 
0.3 < z < 0.6 clusters all with ASCA temperatures (Novicki, Sornig, k, Henry, 
2002). 

Although there is clearly a relation in Figure 1, there is also substantial 

non-statistical scatter. The scatter is approximately Gaussian distributed in 
the logarithm of the luminosity at constant temperature. The dispersion in the 
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Figure 2. The 1 and 2 a contours for a fit of the data in Figure 1 
to L{bol) = C[kT]°'{l + z)^. a. (left) For a cosmology with n^o = 
1.0, J7ao = 0-0. b. (right) For a cosmology with Q^q = 0.3, J^ao = 0.7. 



logarithm of the bolometric luminosity, a"logL(bol), is about 0.2. The scatter 
is reduced if the effects of the cluster cooling centers are removed, either by 
modeling them (Allen &; Fabian, 1998), excising them (Markevitch, 1998), or 
only using clusters that do not have cooling centers (Arnaud &: Evrard, 1999). 

Most observers adopt the simple functional form L{bol) = C{kT]°'{l + Z)^, 
where L{bol) is the bolometric luminosity, to describe their results. The obser- 
vations indicate that a ~ 3. However, simple scaling relations predict a = 2 
and A = 1.5. This a discrepancy is usually resolved by positing that the cluster 
gas has experienced some form of preheating prior to its heating through gravi- 
tational collapse (Evrard Sz Henry, 1991; Kaiser, 1991; Bialek, Evrard, & Mohr 
2001; Tozzi & Norman, 2001; Bower et al., 2001). Preheating raises the tem- 
peratures of low mass clusters more than those of high mass, thereby breaking 
the scaling relations. These preheating models predict little or no evolution of 
the L - T relation, at least to redshifts ^ 1. For example, the model in Evrard 
& Henry (1991) has a = 2.5 and A = 0. An alternative to this idea invokes 
cooling, which is more efficient at lower temperatures, and a small amount of 
post collapse heating to break the scaling (Voit & Bryan, 2001). Little evolution 
is expected in this model as well. Thus the values of a and A provide clues to 
the thermodynamic history of the cluster gas. 

It has recently become possible to constrain the evolution of the L - T 
relation, thereby testing these scenarios. Figure 2 shows the constraints derived 
by Novicki et al. (2002). The evolution is mild; A is consistent with at the la 
level. Similar results have been found by Henry, Jiao, & Gioia (1994), Mushotzky 
& Scharf (1997), Sadat et al. (1998), Donahue et al. (1999), Reichart, Castander 
& Nichol (1999), Fairley et al. (2000), and Arnaud, Aghanian & Neumann 
(2002). Little or no evolution may then imply that the thermodynamic history 
of cluster gas is more complicated than simply falling into the cluster potential 
well, at least for low mass clusters. There was possibly some additional heating 
prior to collapse and/or cooling after collapse. 
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3. Cluster X-ray Luminosity Evolution 

The initial claim for cluster X-ray luminosity evolution was made over a decade 
ago by Gioia et al. (1990) and Henry ct al. (1992) using the Einstein Extended 
Medium Sensitivity Survey (EMSS). The reported evolution was in the sense 
that there are fewer high luminosity clusters in the past. This conclusion has 
not enjoyed universal acceptance. However, there are now many new samples of 
clusters, both at high and low redshift, with which it may be tested. 

3.1. Cluster X-ray Luminosity Samples 

Clusters are not standard luminosity candles. Therefore, searching for luminos- 
ity evolution requires measuring the X-ray himinosity function at two epochs 
and that requires constructing statistically complete samples at the two epochs. 
Fortunately, measuring the X-ray luminosity is relatively easy. It only requires 
~ 25 photons for a 5a luminosity measurement, so large samples are possible 
with relatively modest observing time. 



Table 1. Low Redshift Cluster Luminosity Samples 



Name 


Num 


Sky Cut 


F(0.1,2.4)(cgs) 


ri(dcg^) 


Reference 


BCS 


201 


S > 0.0 


> 4.4 X 10"^^ 


13,578 


Ebeling et al. (1998) 






1 6 |> 20 






REFLEX 


441 


5 < 2.5 


> 3.0 X 10-12 


13,924 


Bohringer et al. (2001) 






1 6 |> 20 









Table 1 gives some properties of low redshift (here defined as z < 0.3) 
statistically complete cluster luminosity samples. Not included in Table 1 are the 
RASSl Bright Sample (De Grandi et al., 1999), which is a subsample of REFLEX 
and the eBCS sample (Ebeling et al., 2000), which adds another 100 clusters at 
F(0.1,2.4) > 2.8 X 10-^2 ergcm-2 s"^ to the BCS but is 75% complete. The 
REFLEX sample was constructed from the second (i.e. fully merged) processing 
of the ROSAT All-Sky Survey (RASS2). 



Table 2. High Redshift Cluster Luminosity Samples 



Name 


Num 


Sky Cut 


F(0.5,2.0)(c 


:gs) 




Reference 


EMSS 


22 


5 > -40 


> 4.6 X 10" 


-14 


778 


Gioia & Luppino (1994) 






1 6 |> 20 










MACS 


119 


^ < 80 


> 6.2 X 10" 


-13 


22,735 


Ebeling et al. (2001) 






S > -40 












1 6 |> 20 










NEP 


19 


a= 18 


> 2.2 X 10" 


-14 


81 


Henry et al. (2001) 






S = 66.5 










160 deg2 


73 


1 6 |> 20 


> 1.3 X 10- 


-14 


158 


Vikhlinin et al. (1998) 


RDCS3 


50 


6 > 20 


> 3.0 X 10" 


-14 


47 


Rosati (2001) 


BSHARC 


12 


6 > 20 


> 1.6 X 10" 


-13 


179 


Romer et al. (2000) 


WARPS II 


78 


6 > 20 


> 6.0 X 10" 


-14 


73 


Jones et al. (2001) 



Table 2 gives some properties of the high redshift cluster luminosity samples. 
Note that most of these samples also have low redshift objects, but only the 
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Figure 3. Ratio of the observed to no-evolution expectation numbers 
of clusters for seven high redshift surveys. The grey bar assumes i^mO = 
1.0, Oao = 0.0; the clear bar i^mO = 0.3, JIao = 0.7. The ratio can not 
yet be made for MACS and WARPS II for the latter case because the 
individual cluster luminosities and redshifts are not yet published. 



numbers oi z > 0.3 objects are listed in the table. The EMSS sample is from 
Einstein Observatory observations, the other samples are from ROSAT. The 
MACS and NEP samples come from the RASS2, i.e. they are contiguous, all 
other samples are from disjoint targeted observations with a region around the 

target excised. 

Typically, low redshift cluster samples contain several hundred X-ray bright 
objects coming from one-third of the sky. High redshift samples usually contain 
several tens of faint objects coming from less than one-percent of the sky. 

3.2. Test For Cluster X-ray Luminosity Evolution In Seven High 
Redshift Samples 

Since the evolution appears to be a lack of objects at high luminosities, simply 
comparing the binned low and high redshift differential luminosity functions will 
not use the full statistical power of the sample. In the limit that the evolution 
is a downward step function above a certain luminosity, such a comparison will 
show no difference between the two epochs (unless upper limits arc plotted). 

A more sensitive test is to compare the observed number of high-redshift, 
high-luminosity clusters to the expected number if the low redshift luminosity 
function did not evolve. This test accounts for clusters that could have been 
found but were not. 
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Figure 4. Ratio of tlic observed to no-evolution expectation numbers 
of clusters for the low redshift components of five high-redshift surveys. 
The sample definition is L(0.5,2.0) > 10"^^ erg s'^ and 0.15 <z< 0.30. 
The grey bar assumes = 1.0, J^ao = 0.0; the clear bar CljnO = 0.3, 
Qao = 0.7. 



Specifically, the test used here compares the number of observed clusters 
with luminosities in the 0.5 - 2.0 keV band > 2 x 10^^ erg s~^ (except MACS, 
whose selection function limits the lowest luminosity to 2.9 x lO^erg s"-*^) and 
with 0.3 < z < 1.0 to the expected number assuming the REFLEX low redshift 
luminosity function did not evolve. Each of the seven high redshift samples 
probes a different region of the L - z plane due to its different flux limit and 
solid angle. The average redshift and luminosity of the clusters after the above 
cuts for all seven samples is 0.495 and 4.7 x 10^^ erg s~-^ respectively (ilmO = 1-0, 
^Ao = 0.0). Figure 3 summarizes the results. All seven samples show evolution, 
most at a statistically significant level. The average ratio of the observed to 
no-evolution numbers is 0.35 ± 0.04 for J^^o = 1-0, ^ao = 0.0, and 0.29 ± 0.04 
for Q^nO = 0.3, Oao = 0.7. 

Figure 3 also shows that one systematic effect, the assumed cosmology, 
is not important. Another systematic effect is the choice of the low-redshift 
luminosity function. Using the BCS instead of REFLEX yields a ratio of 0.38 ± 
0.04 for flrnO = 1-0, J^Ao = 0.0 (the BCS investigators only gave results for this 
cosmology) . Thus the choice of the low-redshift comparison luminosity function 
is also not important. The final systematic investigated was whether the test 
yields a null result, a ratio of 1, at low redshift. Figure 4 shows that it does. 
The average ratio of the observed to no-evolution numbers is 1.03 it 0.08 for 
n^no = 1.0, JIao = 0.0, and 0.87 ± 0.08 for J)„o = 0.3, J)ao = 0.7. This test is 
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Figure 5. X-ray luminosity functions from the low-redshift REFLEX 
(from Bohringer et al., 2002) and combined EMSS, NEP and 160 deg^ 
high-redshift samples. The figure assumes $7^0 = 1-0, r^Ao = 0.0. 



limited to five samples because the selection function for Bright SHARC has not 
been published at low redshift and MACS does not attempt to be complete at 
low redshift. 

These results show that there is no longer a question of whether cluster 
X-ray luminosities evolve, they do. Seven nearly independent surveys show 
some evolution with a combined significance > 15a. Systematic effects do not 
appear to be important. The amount of evolution is about a factor of three at 
L(0.5, 2.0) ~ 4.7 X 10^^ erg s~^ over the redshift range ~ 0.1 to ~ 0.5. It is now 
time to better characterize the evolution. 



3.3. Characterizing Cluster Luminosity Evolution Using Two Large 
High Redshift Samples 

The first large sample comes from combining the 160 deg^, EMSS and NEP 
samples, denoted HEN. There are 109 objects with L(0.5,2.0) > 1 x lO''^ ergs~^ 
and 0.3 < z < 0.85. This size is comparable to the low-redshift samples. There 
are no clusters in common among the three samples and the overlap of their 
survey regions on the sky is small, only ~ 5%, which has been accounted for 
statistically. The HEN sample is large enough that a simple comparison of the 
luminosity functions from it and REFLEX demonstrates the evolution clearly, as 
Figure 5 shows. A Schechter function fits both the REFLEX and HEN samples 
well. Figure 6 gives the 1, 2 and 3 a contours for the slope and characteristic 
luminosity, again demonstrating evolution at high statistical confidence. At least 
the Schechter function normalization and slope are evolving, but the charcteristic 
luminosity may not be. Comparison of HEN with REFLEX shows a factor of 2 
evolution at L(0.5,2.0) = 2 x 10^^ ergs~^ over the redshift interval 0.11 to 0.45 

{nmo = 1.0, = 0.0). 
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Figure 6. The 1, 2 and 3a contours for the slope and characteristic 
luminosity of individual Schechter function fits to the REFLEX and 
HEN samples, a. (left) For a cosmology with = 1-0, J^ao = 0.0. b. 
(right) For a cosmology with = 0-3, ^ao = 0.7. 

It is the nature of flux-limited samples that different luminosity objects are 
predominately found at different redshifts. Thus the luminosity functions in Fig- 
ure 5 do not refer to unique epochs. This situation may be overcome to some ex- 
tent by fitting to an evolving luminosity function, termed the AC model. The AC 
model generalizes the Schechter funtion to: n(L, z) = no(z)L~"(^)exp(— L/L*), 
with no(z) = no(l -|- z)^ and a(z) = ao(l + z)^. A maximum likelihood fit of the 
AC model to the unbinned luminosity and redshift pairs for the REFLEX and 
HEN clusters gives the following: no = 4.651^;^^ x 10-'^Mpc-^(10^^erg s-i)°-\ 
L* = 6.86lti X lO^^crg s'S ao = -(1.70t{|:if ), A = -(LIT ± 0.84), C = 0.90t[|;^ 
(r^mO = 1-0, i^Ao = 0.0). The errors are at 68% confidence for 4 parameters with 
the error on uq determined from letting the errors on the other 4 parameters take 
their extreme values. Figure 7 gives the best fitting AC luminosity functions, 
now at two particular epochs. 

The second large high redshift sample is the Massive Cluster Survey or 
MACS. This survey seeks a large sample of high luminosity clusters at high 
redshift. About 300 objects are expected if there is no evolution. Presently the 
survey is ~ 90% complete and has 119 clusters. Note that in Section 3.2 an 
additional 10% systematic error has been added in quadrqture to the MACS 
statistical errors. 

A final sample of ~ 130 objects means that MACS detects evolution also 
at high significance. Figure 8 shows a nonparametric representation of this 
evolution, comparing the MACS and eBCS integral luminosity functions. The 
comparisons in Section 3.2 folded a model fit to the REFLEX luminosity func- 
tion through each survey selection function. MACS finds a factor of ^ 10 
evolution at L(0. 1,2.4) > 10^"'' crgs^^ over the redshift interval 0.05 to 0.55 
{^rnO = 1-0, f^AO = 0.0). Figure 8 also shows that the AC model fits the MACS 
data well, with no adjustable parameters and with few clusters in common be- 
tween eBCS/MACS and REFLEX/HEN. 

Two high statistics, high redshift samples agree on the amount and nature 
of cluster X-ray luminosity evolution. The previous apparent disagreements were 
probably just poor statistics. 
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Figure 8. A nonparametric description of X-ray cluster luminosity 
evolution detected in the eBCS/MACS samples. The integral luminos- 
ity function is shown as a function of redshift. Strictly speaking, the 
MACS results are lower limits, but the estimated completeness of this 
survey is ~ 90%. Also shown is the AC model with its la error band. 
There arc no adjustable parameters used to fit the AC model to the 
eBCS/MACS. The figure assumes nrnO = 1-0, Oao = 0.0. 
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3.4. Constraining Cosmology From Cluster Luminosity Evolution 

Given the discussion in the Introduction, the question naturally arises whether 
cluster luminosity evolution may be used to determine the cosmology of the uni- 
verse. The difficulty is the extra evolution experienced by the gas in addition 
to that driven by gravity, as the L - T relation shows. Since X-ray luminosity is 
proportional to the square of the gas density, will this extra evolution be large 
enough to mask that from gravity? The weak evolution of the L - T relation 
shows that the extra evolution occured before what is currently observable, pos- 
sibly allowing a simple correction. Comparing data to a theory that incorporates 
the observed L - T (or L - M) relation into the gravitational evolution predic- 
tions has in fact been used to constrain cosmological parameters from cluster 
luminosities (eg Henry et al., 1992; Borgani et al., 2001; Viana, Nichol &; Liddle, 
2002). 

4. Cluster X-ray Temperature Evolution 

Perhaps a cleaner way to measure the evolution of cosmic structure using X-ray 
observations of clusters is via cluster temperature evolution. The temperature 
is more likely to reflect better the underlying gravity driven evolution than the 
luminosity. 

4.1. Cluster X-ray Temperature Samples 

Clusters are not standard temperature temperature baths. Just as with luminos- 
ity evolution, measuring temperature evolution requires constructing tempera- 
ture functions at two epochs. Again, statistically complete samples are required 
at the two epochs. However it is substantially more difficult to measure cluster 
temperatures than luminosities. At least 1000 photons are required to provide 
a reasonable temperature measurement. Obtaining large samples is difficult. 
Obtaining even two independent samples at high redshift has so far not been 
possible. Table 3 gives some properties of cluster temperature samples. The last 
column of the table gives the numbers of clusters with measured temperatures 
from ASCA and with temperatures estimated from the L - T relation. 



Table 3. Cluster Temperature Samples (z<0.3 All Sky; z>0.3 778 deg^) 



Reference 


Num 


z Cut 


T Cut 


F(0.1,2.4)(c 


gs) 


ASCA/L-T 


Edge et al. (1990) 


46 






> 1.7 X 10- 


■ 11 


0/2 


Henry & Arnaud (1991) 


25 


< 0.17 




> 3.0 X 10- 


11 


25/0 


Markevitch (1998) 


30 


0.04-0.09 




> 2.0 X 10- 


11 


26/4 


Blanchard et al. (2000) 


50 


< 0.10 




> 2.2 X 10- 


■ 11 


34/0 


Pierpaoli et al. (2001) 


38 


0.03-0.09 


> 3.5 


L-T, z 




?/l 


Ikebe et al. (2002) 


61 




> 1.4 


> 2.0 X 10- 


■ 11 


56/2 


Henry (2002) EMSS 


22 


0.30-0.85 




> 1.2 X lo- 


13 


21/1 



The low redshift samples are comprised of bright well-known clusters and 
are all more or less the same. The Edge et al. and Henry Sz Arnaud samples 
come from non-imaging surveys, the rest of the low redshift samples come from 
ROSAT. The Markevitch and Ikebe et al. samples have corrections for cooling 
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Figure 9. A nonparameteric, unbinned presentation of the cluster 
integral temperature function from the Henry & Arnaud (1991) and 
EMSS samples (Henry, 2002). The figure assumes = 1-0, r^Ao = 
0.0. 



centers; about 20% of the Blanchard et al. sample has these corrections. The 
high redshift clusters are unresolved by ASCA, which, when combined with the 
low statistics, does not permit corrections for cooling centers. Figure 9 shows 
the integral temperature functions at two epochs. The same factor ~ 2 evolution 
seen in the luminosities is also seen in the temperatures. 



4.2. Cosmological Results From Cluster X-ray Temperature Evolu- 
tion 

Deriving cosmological constraints from X-ray observations is a two-step process. 
The fundamental variable of the theory is an object's mass so the process starts 
with a mass function. The analytic mass function that best matches the large 
N-body simulations comes from Sheth and Torman (1999). This function agrees 
better with the simulations than does the venerable Press - Schechter (1974) 
function. However, the differences are only a factor of ~ 2. 

The evolution of the mass function is very sensitive to cosmological pa- 
rameters, but the mass of a cluster is difficult to determine. So, the mass 
must be converted to a more easily observed quantity, in this case the tem- 
perature. A top hat collapse is most often used to do that. This theory 
gives kT = 1.42/3TM[f^moA(Omo, ^Ao, z)] V3(hMi5)2/3(l + z). Here /?tm is a fac- 
tor near unity that accounts for the effects of partial virialization of the cluster. 
Its exact value comes from liydrodynamical simulations. A(r2mOi ^AO) z) is the 
ratio of the average cluster mass density to the background mass density for a 
cluster that collapses at redshift z. It has recently been determined for arbitrary 
^mo and Oao by Pierpaoli et al. (2001). M15 is the cluster mass in units of 
1O^^M0. 
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Figure 10. The 1 and 2 a contours for and Oao from cluster tem- 
perature evolution, cosmic microwave background and supernovae. The 
supernovae and CMB constraints are from Perlmutter et al. (1999a) 
and Efstathiou et al. (2002) respectively. 



One final issue concerns the selection function. X-ray clusters are selected 
by flux, not temperature. That is the solid angle surveyed to a given flux limit 
is known. Pierpaoli ct al. (2001) circumvent this by using the L - T relation 
to insure that their temperature and redshift selection was compatible with the 
flux limits of existing samples. Usually, the flux limit is recast as a limit in the 
luminosity-redshift plane and then the luminosity is converted to a temperature 
via the L - T relation. The dispersion in the L - T relation is included by 
averaging over it which gives the solid angle surveyed for a given cluster of 
temperature kT at redshift z: 

^(kT,z) = / ^^#fS%T«-P - (log(C[kT]") - logL)V2a2(logL) 
J -\/27rcr^(logL) 

There has been a great deal of recent work in this area. Among the many 
results are: BahcaU & Fan (1998), Donahue & Voit (1999), Henry (2000; 2002), 
Pierpaoh et al. (2001), and Viana & Liddle (1999). 

Constraints on cosmological parameters presented here come from a max- 
imum liklihood fit of the evolution theory to the unbinned temperature and 
redshift data for the samples in Figure 9. There are four interesting parameters 
in the theory: the matter and dark energy densities and the shape and nor- 
malization of the matter power spectrum. The effects of temperature errors are 
included in the likelihood function. Figures 10 and 11 show some of the results 
(Henry, 2002). 

Figure 10 shows that clusters provide a band of constraints in the 0,„io — JIao 
plane of comparable width but different orientation to those from supernovae 
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Figure 11. a. (left) Constraints on cjg and OmO shown as 1 and 2 
a contours from cluster evolution and a 2 a contour from the CMB 
(Melchiorri & Silk, 2002). The CMB constraint assumes h = 0.71 ± 
0.07 and $7^0 + ^Ao = 1. b. (right) Constraints on w and ilmo shown 
as 1 and 2 a contours from cluster evolution and supernovae and a 1 cr 
contour from the CMB. The supernovae and CMB constraints are from 
Perlmutter, Turner k White (1999b) and Bean k Melchiorri (2001), 
respectively. All results assume that QmO + ^Qo = 1 and the CMB 
result further assumes h = 0.72 ± 0.08. 

and microwave background observations. Further, all three bands intersect at 
the same point in the plane. 0,mo = 0.4, IIao = 0.6 is consistent with all three 
data sets at the ~ la level. For example, Rubino-Martin ct al. (2002) show 
that all CMB data reported to date, plus the supernovae data from Perlmutter 
et al. (1999a), plus assuming h = 0.4 - 0.9 yields flmO = 0.35 ± 0.04 and 
Qao = 0.68 ± 0.04 at 68% confidence. 

The normalization of the power spectrum, dg, agrees with the microwave 
background normalization evolved to the present. Figure 11a. Thus the fluctu- 
ations seen in the CMB really do grow to form the clusters seen today. The 
present cluster sample size is not large enough to provide tight constraints on w. 
Figure lib. However, w = -1 is a good bet considering the constraints provided 
by supernovae and the CMB. 

4.3. Cluster X-ray Temperature Evolution Systematics 

It is of course important to consider the effects of systematic errors on the 
constraints derived in the last section. The three effects considered here are: 1. 
Using the older Press - Schechter (1974) mass function; 2. Using an empirical 
normalization of the mass - temperature relation, i.e. /3tm, from Finoguenov, 
Reiprich &; Bohringer (2001) instead of a theoretical normalization; 3. Setting all 
temperature errors to 2% or the actual value, whichever is less. All systematics 
effects considered change the constraint on by ~ la, the constraint on ilmo 
by ~ 0.5(7 and the constraints on JIao and w by <IC la. Thus statistical errors 
appear to dominate the error budget with the present sample size of about fifty 
objects. 
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5. Conclusions 

The L - T relation shows that the hot X-ray cluster gas experienced a more com- 
plicated thermodynamic history than that resulting from the cluster collapse. 
There was possibly some heating prior to collapse and/or cooling with lesser 
heating after collapse. 

There is no longer a question whether cluster X-ray luminosities evolve, 
they do. Seven nearly independent surveys show some evolution and system- 
atics do not appear to be a big effect. Large high redshift samples, i.e. those 
containing more than 100 clusters, are needed to characterize the evolution. Two 
such samples currently exist. HEN shows a factor of 2 evolution compared to 
REFLEX at L(0.5,2.0) = 2 x lO''^ ergs"^ over the redshift interval 0.11 to 0.45 
(S7m0 = 1-0, f^Ao = 0.0). MACS shows a factor of 10 evolution compared to cBCS 
at L(0. 1,2.4) > 10"^^ ergs-i over the redshift interval 0.05 to 0.55 {nmo = 1.0, 
Qao = 0.0). 

A cosmology with = 0-4, ^ao = 0.6 is consistent with cluster X-ray 

temperature evolution, the cosmic microwave background and supernovae at 
the ~ la level. The systematics of cluster temperature evolution constraints 
on cosmology are a minor concern at present, i.e. they are about the same size 
as the statistical errors. However, this method is so statistically powerful that 
systematics will be the dominate effect as soon as the sample sizes become only 
a factor of two or three larger, at least for the determination of ag- 

A more positive statement would be that an X-ray cluster survey of 10,000 
deg2 to a flux of F(0.5, 2.0) > 5 x 10"^'^ erg cm'^ s"^ would yield ~ 18, 000 clus- 
ters to z ~ 1.5 and provide constraints of similar statistical quality as the up- 
coming Planck and the proposed SNAP missions (Pctrc ct al., 2001). In fact 
the great complementarity of cluster evolution, cosmic microwave background, 
and supernovae constraints exhibited in Figure 10 would provide a check on the 
systematics of all three methods. 
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